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Abstract. We discuss different approaches to the problem of reproducing the 
observed features of nuclear single-particle (s.p.) spectra. In particular, we analyze 
the dominant energy peaks, and the single-particle strength fragmentation, using the 
example of neutron states in 208 Pb. Our main emphasis is the interpretation of that 
fragmentation as due to particle- vibration coupling (PVC). We compare with recent 
Energy Density Functional (EDF) approaches, and try to present a critical perspective. 



1. Introduction 

One of the most striking features of nuclei is the validity of the shell model where the 
nucleons move independently in a smooth single-particle (s.p.) potential. However, self- 
consistent mean field (SCMF) approaches [1] are much more sucessful in reproducing the 
bulk properties that emerge from the mean field picture (like masses, that can nowadays 
be reproduced by Skyrme-HFB calculations with a remarkable accuracy of about 0.6 
MeV [2]) than in accounting for the s.p. states. 

If SCMF is viewed as an approximate realization of an exact, yet unknown and 
perhaps unfeasible, Density Functional Theory (DFT), it is not surprising that quantities 
like the total energy are well accounted for since they belong to the domain of this theory 
based on the well-known Kohn-Sham theorem. (We do not enter here the problem if 
this theorem can be expected to hold in a self-bound system since this discussion can 
be found elsewhere in the present volume). 

If single-particle states are not, strictly speaking, within the DFT framework, two 
possible routes can be undertaken. Some authors are currently trying to improve 
the accuracy of present DFT implementations aiming at functionals with so-called 
"spectroscopic" accuracy (referring to, of course, single-particle spectroscopy). Another 
point of view, which is the preferred one by the authors of the present contribution, 
is to generalize the shell model to the so-called "dynamical" shell model. It has to 
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be noted that the former attitude does not solve the problem of evaluating the so- 
called spectroscopic factors or, more generally, the well-known fragmentation of the s.p. 
strength. 

The potential emerging from SCMF calculations is static, but in general nonlocal in 
space (or, in equivalent words, velocity- or momentum-dependent). The fluctuations of 
this average potential lead to collective modes, in particular surface vibrations. Taking 
into account the coupling of these modes to the s.p. motion, the shell model acquires 
a dynamical content. Thus, the average potential becomes also nonlocal in time, being 
characterized by an energy (or frequency) dependence. 

The velocity dependence of the s.p. potential can be characterized by the so-called 
/c-mass rrik, and its frequency dependence by the w-mass m w . Most experiments only 
probe the product of these two quantities, that is, the "total" effective mass m* (see 
the definitions in Sec. [2]). The dynamics of the shell model affects different nuclear 
properties as: the fragmentation and related spectroscopic factors of the s.p. states, 
their density (which is proportional to m* near the Fermi energy), the s.p. spreading 
widths, and the imaginary component of the optical potential. A unified description of 
the s.p. motion at positive and negative energy eventually emerges. 

The importance of the dynamical couplings was first demonstrated in the case of 
electrons in metals and of normal liquid 3 He. These results inspired the work of Bertsch 
and Kuo [3] on the enhancement of the effective mass near the Fermi energy in finite 
nuclei. The review article [I] contains the formalism reported in the next Section as well 
as the many results obtained in the 80s. The expression "dynamics of the shell model" 
is after G.E. Brown [5]. 

The old calculations of the 80s are mostly not consistent. It is hard to extract 
quantitative conclusions because of the many approximations involved. More consistent 
calculations have been recently done. They can be confronted with modern EDF-based 
approaches. We shall do it in the simple case of 208 Pb, showing that the problem of the 
s.p. states is still an open one. 

2. Formalism of the dynamical shell model 

One can introduce using standard many-body techniques a "mass operator" J\4, defined 



by 



[u + {2m)- l V 2 r 



'] G(r, r"; u) = 8(f- f) + / d 3 r' M(r, r*\ u)G(r, f*; u) (1) 



in terms of the one-body Green function (Lehmann representation) 





Here 



Mrl = (4 A - 1) \a^\4 A) ) 1 
Mrl = (4 A) H^ i p A+1) }- 



(3) 
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ifrif is a suitably normalized eigenstate of the Hamiltonian of the (A-l)-particle system 



with energy E- 



(A-l) j(A+l) . 



li 



is a suitably normalized eigenstate of the Hamiltonian of the 



(A+l)-particle system with energy and 
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Eqs. (CED and (ED give 

[w + (2m)' 1 V 2 ] <j) p {r) 



d 3 r' Ai(f, r'; to, 



pJYp 



0. 



(4) 



(5) 



This s.p. wave equation has two important properties: i) It has discrete eigenvalues 
corresponding to the energies of the bound states of the (A+l) and (A-l) systems, ii) 
In the continuum gives the exact diagonal (elastic scattering) element of the S matrix 

The mass operator, in what follows, will contain the Hartree-Fock (HF) mean field 
and the contribution of the coupling of the s.p. HF states to the vibrational modes 
(particle- vibration coupling or PVC), in particular to collective modes obtained from 
self-consistent linear response theory (Random Phase Approximation or RPA). 

Although one could solve (JH]) directly, the usual way has been so far to solve it 
at the HF level and add the contribution of the particle-vibration part of the mass 
operator within second-order perturbation theory. The expressions for this second-order 
correction can be found in most of the references quoted herein. 

To work out the expressions for mf., and m*, it is easier to consider a uniform 
system [3]. Thus, having the relation 



(6) 



between the energy E and the momentum k of the quasiparticle (dressed s.p.), and V 
being the real part of the complex mass operator, the effective mass is defined 

dE k ^ 



dk m* 
Then, it is easy to write 



to obtain 



m 
m 

m* 
m 



1 + 



m 



dV{k-uj) 
dk 



uj=E(k) 



u=E(k) 



m 



m 



(8) 



(9) 



3. The contradicting results 

3.1. Particle-vibration coupling: old and new results 

In Tables 4.3a and 4.3b of Ref. [3], an extensive review of the results obtained for m w in 
208 Pb by nine groups in the period 1968-1983 can be found. Rather different frameworks 
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had been adopted, no one being fully self-consistent: s.p. potentials range from harmonic 
oscillator (HO) to Woods-Saxon (WS) or HF with Skyrme forces; residual interactions at 
the particle-vibration vertex are either multipole-multipole forces, or forces of Landau- 
Migdal type, or Skyrme forces but with velocity- dependent terms dropped, or even 
G-matrix interactions. 

The average value for is in the range 1.2-1.4 for the neutron valence shells, 
and in the range 1.3-1.6 for the proton case. The inverse of gives the spectroscopic 
factor reported in table 4.4 of Ref. [I]. In general, it has been impossible to find an 
accurate agreement between theory and experiment, either at the level of main energy 
centroids of the s.p. strength distribution, or at the level of spectroscopic factors. The 
problem of spectroscopic factors had been not considered too seriously, in view of the 
large experimental errors and ambiguities associated with them. Nowadays, new analysis 
are being carried out [7], but the theoretical ambiguities associated with the extraction 
of spectroscopic factors are still under discussion (see e.g. the contribution by R. J. 
Furnstahl and A. Schwenck in the present volume). 

If we concentrate on s.p. energies, the most consistent among the calculations 
reported in [I], namely that by V. Bernard and N. Van Giai [8], still underestimates the 
density of levels. 

Table 1. Shifts (that, is differences between dressed and bare single-particle energies) 
associated with the dominant s.p. neutron components in 208 Pb. See the text for a 
discussion of how they have been obtained. 
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The most recent calculations for 208 Pb in the framework of dynamical shell model, 
or PVC model, have been reported in Refs. [9] and [ID] . The results are displayed in 
Table [I] and [2] and compared with the results of Ref. [11] (cf. below) and with the 
experimental energies, in the case of the neutron valence particle and hole states. In 
Table[T] we display, for Refs. [9] and [10] . the differences between the s.p. energies £j 
obtained with PVC included, and the mean field energies ef\ These differences A£j 
will be called "shifts" in what follows. (See below for a definition of the shifts from Ref. 
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in])- 

In the covariant theory of Ref. [9], the mean field is associated with the NL3 
[12] effective Langrangian, which is also employed to obtain the phonons within the 
relativistic RPA (RRPA). The Dyson equation ([2]) for the s.p. Green's functions is 
solved in the diagonal approximation for the mass operator Ai. A rather simple particle- 
vibration model is used, in which the coupling vertex is obtained with the matrix element 
of the residual interaction times the transition density of the vibrations. 

The calculations of Ref. [10] are performed in the framework of Skyrme HF, with the 
PVC contributions to the mass operator Ai calculated within second-order perturbation 
theory. The phonons are obtained through fully self-consistent RPA, and all the terms 
of the particle-hole (p-h) interaction are kept also in the PVC vertex. There are strong 
cancellations due to the velocity dependent terms that lead to small shifts, especially 
for the hole states. The calculations have been done using the SLy5 Skyrme set [13J. 
Calculations with other Skyrme functionals are in progress. 

Table 2. Dressed single-particle energies and spectroscopic factors S associated with 
the dominant s.p. neutron components in 208 Pb, compared with the experimental 
findings. 
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In Ref. [11] . the HF problem is solved for the even-even and even-odd (A±l) 
systems using the SLj4 L interaction. The "dressed" energies are obtained from an 
equation like the present Eq. (jlj). In this case Eq is the total energy of the even-even 
system and E p (Eh) is the total energy of the odd system when a given occupancy 
(vacancy) is imposed. The "bare" energies are obtained by excluding the so-called 
mass, shape and spin polarization effects: in practice this corresponds, respectively, 
to neglecting the proper change of the center-of-mass correction, to imposing spherical 
symmetry, and to discarding time-odd terms. 

One can see from Tabled] that shifts obtained by RMF plus PVC and Skyrme plus 
PVC are different, and it is hard to connect either of them to the calculations of [TT]. If 
we look now at the final "dressed" energies obtained in the different approaches (Table 
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[2]), and compare with experimental values, we conclude that the open problem is that 
even for the "benchmark" nucleus 208 Pb we do not dispose of a consistent calculation 
which reproduces well the experimental data. This is also true for spectroscopic factors. 

In the calculations of the spectroscopic factors, the effects of the short-range 
correlations (SRC) should also be included. On the one side, in a recent work [H], it 
is confirmed that the PVC involving long-range correlations (LRC) can be identified as 
the main mechanism, the SRC explaning only a small fraction (up to about 10%) of the 
deviations of the spectroscopic factors from the independent s.p. values. Consistently, 
from a qualitative point of view, in Ref. [I] it was concluded that the main contribution 
of highly excited core polarized states to m u , and therefore to the spectroscopic factors, 
can be considered as a constant in the vicinity of the Fermi energy. The subtracted 
dispersion relation technique, with different choices for the behaviour of the imaginary 
part of the mass operator at high energy, was used. 

3.2. EDF extensions 

The lack of an accurate description of s.p. states is also clear for EDF practitioners 
not only in keeping with the findings of Ref. [11]. The authors of [15] have concluded 
that standard forms of the Skyrme energy density functional do not allow that accurate 
description. The effect of tensor forces has been studied very carefully by the authors of 
|16[ ITT], but their conclusion is that the tensor terms, albeit important, cannot remedy 
the deficiencies of the central terms: in practice, they cannot produce a functional which 
is more accurate as far as s.p. states are concerned. 

In Ref. [H], gradient corrections up to next-to-next-to- next-to-leading order 
(N 3 LO, in practice derivatives up to sixth power) have been added to a local EDF, 
in order to better reproduce surface effects, for example on the effective mass m w (cf. 
also [19]). The precision on s.p. states is improving [20] . 

4. Conclusions 

We still lack, despite the tremendous progress in many aspects of nuclear structure 
physics, an accurate description of s.p. strength distributions. While new radioactive 
beam facilities have already started to show how this kind of observables evolve with the 
decrease of the neutron or proton separation energies [2TJ [22] , even for standard nuclei 
we have troubles to get accuracy for dominant s.p. energies and spectroscopic factors 
(using models that have had big success in reproducing other kind of observables). 

The addition of higher power gradient corrections to a local EDF may improve 
the s.p. energies. We deem, however, that a dynamical theory is necessary to describe 
quantities like the spectroscopic factors. The merging of the two approaches defines 
anyway, to some extent, the future of this fundamental field of nuclear structure research. 
In a many-body approach like the dynamical shell model one believes that the basic 
object should be a time- or energy-dependent quantity like the s.p. Green's function. 
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EDF practitioners aiming at functionals with spectroscopic accuracy believe that a 
static approximation, in which the Green's function reduces to a static density, may 
be a sound approximation. This issue asks for a clear answer based on comparision of 
results without uncontrolled approximations. 
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